Integration Theory for infinite dimensional 
volatility modulated Volterra processes 

Fred Espen Benth*and Andre SiiB^ 
March 29, 2013 



We treat a stochastic integration theory for a class of Hilbert-valued, 
volatihty-modulated, conditionahy Gaussian Volterra processes. We apply 
techniques from Malliavin Calculus to define this stochastic integration as a 
sum of a Skorohod integral, where the integrand is obtained by applying an 
operator to the original integrand, and a correction term involving the Malli- 
avin derivative of the same altered integrand, integrated against the Lebesgue 
measure. The resulting integral satsifies many of the expected properties of 
a stochastic integral, including an Ito formula. Moreover, we derive an al- 
ternative definition using a random-field approach and relate both concepts. 
We present examples related to fundamental solutions to partial differential 
equations. 

Keywords: Stochastic integration, Volterra processes, Gaussian random fields, Malli- 
avin Calculus 

1 Introduction 

Let throughout this article < T < oo be a finite time horizon, fix t G [0, T] and let 
Til, 7^2 and be three separable Hilbert spaces. As the main object of investigation 
in this article we introduce the following process 

X{t) = I g{t,s)a{s)6B{s), (1.1) 
Jo 

where i? is a cylindrical Wiener process on Hi, a is stochastic process on a time interval 
[0, T] with values in L{'Hi,'H2), not necessarily adapted to the Wiener process B, and g is 
a deterministic function depending on two time parameters such that g{t, s) G ^(^2i ^2) 
for all < s < t < r. In order for the stochastic integral to be well-defined one has to 
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assume that g{t,-)a{-) is Skorohod integrable on [0,t] so that X(t) exists as a random 
element in 7^2- The aim of this article is to define a stochastic integral with respect to 



where we assume that Y{t) € L{T-i2,T~iz) for all t € [0,T]. With the integration concept 
we develop, we will see that naturally Z{t) € 7^3 for all t E (0,T]. 

We want to point out some remarkable facts about the objects we have just introduced. 
First, g{s,s) does not have to be defined for any s G [0,T], it can be singular at these 
times. Note moreover that one could put suitable measurability conditions on a, such 
as predictability, but in general this is not necessary. In such a case the integral in (jl.ip 
would turn out to be an Ito integral in a Hilbert space and the condition for the existence 
of the integral would be 



At this point we note that our considerations go beyond the classical semimartingale case 
as treated in |17] for real-valued and [7] for Hilbert- valued (semi)martingales. We will 
see however in Section [4.21 that under some conditions on g (in particular that g{s,s) 
exists for all s € [0,T]), X will turn out to be a semimartingale and that the integral 
in ()1.2p and the classical integral with respect to a semimartingale, denoted hy Y ■ X 
coincide in some cases. In fact, if X is a semimartingale, then the difference between 
these two integrals can be compared to the situation of the Ito and Skorohod integral 
with respect to the one-dimensional Brownian motion. Next we list some examples in 
order to show the wide range of processes X that can be used as integrators in (|1.2|) . 

Example 1.1 (Ambit fields). The situation that motivates our problem of defining (jl.2p 
comes from ambit processes, see 0/. There we deal with random fields consisting of a 
stochastic integral over a random field defined as follows: 



where D{t,x) C , W is a Gaussian noise white in time and (possibly) correlated in 
space, and g and a are deterministic and random fields, respectively. These processes are 
included in our setting if interpret these equations in the Hilbert space sense of /77 where 
7i2 is interpreted as Lp'{D{t,x)) and X(t, •) is assumed to be in some L? with respect 
to the spatial parameter. In Section we come back to this example and also derive an 
integral with respect to a random field, similar to 119]. 

Example 1.2 (Gaussian Processes and VMBV). A one- dimensional subclass of X has 
already been treated in and more generally in In the former paper, the au- 
thors study integration with respect to M-valued Gaussian processes where they assumed 



the stochastic process X = (X(t))tG[o,T] i that is we want to derive an integration theory 
for the integral 





ds < oo. 
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a{s) = 1 for all s € [0,T]. In the latter paper the authors considered the possibility 
of a non-trivial a and they referred to those processes as in (jl.ip in one dimension as 
volatility modulated Volterra processes driven by Gaussian noise (VMBV). In this article 
we generalize both Gaussian and VMBV processes to infinite dimensions. A particular 
example is fractional Brownian motion in infinite dimensions. Ghoose for all s G [0,T], 
a{s) = Q^l"^ where Q is a nonnegative, selfadjoint, trace-class operator, let H G (0, 1) 
and set 

g{t,s) = CH{t - sf-'/' + ^) - sf-'/'{l - {s/u)'/'-^)du. 

Then X is a Hilbert-valued fractional Brownian motion with Hurst parameter H and 
covariance operator Q. 

Example 1.3 (Solutions to S(P)DE). Another application is stochastic integration with 
respect to the solution of a stochastic differential equation in a Hilbert space. This in- 
cludes solutions to SPDEs interpreted in the sense of J7l. Let for instance X be the mild 
solution of 

dX{t) = -AX{t) +a{X{t))dB{t), 

with Xq = 0, where A is an unbounded linear operator, a is a deterministic function 
subject to some regularity conditions and B is again a cylindrical Wiener process on 
some Hilbert space, see |^ Chapter 6] for a detailed treatment of these equations. Then 
if —A generates a strongly continuous semigroup of linear operators (fi'(t))te[o.r]; mild 
solution to this equation is given by the following integral equation 

X{t)= f g{t-s)a{X{s))dB{s), 
Jo 

which has the form of (jl.ip . With the help of the theory we develop in this article, we 
are then able to define a stochastic integral with respect to this solution X . A particular 
example here are Ornstein- Uhlenbeck processes in infinite dimensions given by the SDE 
dX{t) = —AX{t) + FdB{t), where F is a bounded linear operator. Let the Co-semigroup 
generated by —A be denoted by g{t, s) = exp(— (t — s)A) so that in this case 

X{t)= [ exp{-{t- s)A)FdB{s). 
Jo 

The paper is structured in the following way. In Section [5] we list some fundamental 
results and in the subsequent Section [3] we give the motivation for the definition of our 
integral. Sections 14.11 and [S] are dedicated to showing properties of the integral and an 
Ito formula. In the final Sections [6] and [7] we provide a random-field integration approach 
to the integral and show their equivalence. 

Throughout this article C denotes a positive generic constant, which may change from 
line to line without further notice. 
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2 Prerequisites 



2.1 Vector Measures 

This section deals with the generahzation of measures to set functions taking values in 
a Banach space, so called vector measures. The case which is most important for us 
is when the vector measure is defined on subsets of M+ taking values in the space of 
the bounded linear operators L{'H2,7i2), where 7^2 is a separable Hilbert space. After 
providing the definition and the important concept of total variation of a vector measure, 
we list the most relevant properties here and refer to [S], [5] or [T2] for more details. 

Definition 2.1. Let (J^,,^) be a measurable space and let B be a Banach space. A 
set function fx : J- B is called a finitely additive vector measure, or in short vec- 
tor measure, if U F2) = A'(-^i) + /^(-F2) for any two disjoint sets Fi,F2 € 
Moreover if, for any sequence {Fn)nm ^ of pairwise disjoint subsets of F, we have 
/i(|J^]^F„) = Yl'^=i f^i^n) , then /i is called a countably additive vector measure. Note 
that the convergence of the sum takes place in the norm topology of B. 

The (total) variation of a vector measure is the set function on (F,^) with values 
in M+ U {00} defined for all F ^ ^ by 

WF) :=sup 

where the supremum is taken over all partitions it of F into a finite number of pairwise 
disjoint sets A £ F. If |^|(^) < 00, then ^ is said to be a vector measure of finite 
variation. 

The total variation of a vector measure is the smallest of all nonnegative, additive set 
functions A such that < \{F) for all F G ^. If a countably additive vector 

measure has finite variation, then is also countably additive. For vector measures (not 
necessarily having finite variation) an integration theory similar to that for M+-valued 
measures can be developed, see [9, Section III]. 

We briefly list some properties of the Lebesgue-Stieltjes integral with respect to a 
Banach-valued function which are important in the remaining paper. Given a ;B-valued 
function 5 on a finite or infinite interval of R, which is assumed to have bounded variation. 
Then we can define a vector measure on all finite subintervals [a,b] by //^([a, 6]) := 
g(b) — g{a). Then we extend the measure //^ onto the Borel cr-field of M+ as its Lebesgue 
extension. We denote the extension by fig, too. Next, we fix a function / : — )• 
L{B,Bi), where Bi is another Banach space. Then, if / is a yU^-integrable function over 
an interval [a,b], then we denote the integral with respect to fj,g throughout this article 
by f{s)g{ds) instead of J^^ f{s)ng{ds). A special choice for g is the identity on M. In 
this case the integral is also known as Bochner integral or Pettis integral depending on 
measurability properties of / and integrability properties of ||/||, see [51 Chapter II]. 

Finally, we provide a notion of absolute continuity of the vector measure fig with 
respect to the one-dimensional Lebesgue measure A. For this, consider the measure 
space ([0, T], =^([0, T]), A). Let g : [0,T] B he a function of finite variation and 
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assume that there exists some function 4> : [0,T] — > i3 such that g{t) = <f>{s)ds. This 
holds for instance when g is Frechet differentiable with derivative (/). Assume moreover 
that / is a /^g-integrable function. Then we have that for every A E ;S([0,T]) 

/ f{s)gids) = [ /(s)/u,(ds) = / f{s)cPis)ds, 
J A J A J A 

see [21 Section III. 11] for more details. Throughout the rest of this article we will apply 
the facts stated in this section to the measure generated by the L('H2, 'H2)-valued function 
g used in (jl.ip . always integrating with respect to the first time argument while leaving 
the second one fixed. 



2.2 Multidimensional Stieltjes Integration 

In this section we give a quick reminder about an extension of the integration theory 
treated in the previous section. We want to define the Lebesgue-Stieltjes integral with 
respect to a function that has more than one argument, see [9l Chapter VII]. For this 
we need the concept of (locally) bounded variation in the case of functions with several 
variables, so called BV functions. Let for this [/ C M'^ be an open subset and let Oc{U) 
be the set of all precompact open subsets of U . Then a function g : — )• M is said to 
be of locally bounded variation in U if g £ Lj^ci^) ^^^'^ 

V{g,0) :=sup| / g{x)^(t>{x)dx;4>eCl{0;R''),Uh^^u)<^ 

is finite for all O € Oc{U). Then the set of all functions with (locally) bounded variation 
forms a Banach space, which is nonseparable. For our needs, the most important prop- 
erty of BV functions is that they are precisely those functions with respect to which one 
can define a Stieltjes integral of all continuous functions. So one defines an integral like 
that by starting with the simple functions f{x) = lAix) = Y['j=i l[aj,b,] (^^i)) i-^- indicator 
functions of sets A = X'^^i[aj,bj], where a, 6 € R*^ and a < b (coordinatewise). For 
these functions we define the Stieltjes integral with respect to g and define the notation 
9iA) by 



9{A) :- 



[ f{x)g{dx) = j2i-^y E 9{x). (2.1) 

JM.d fr'a 



Xie{bi,ai};i=l,... 
\{i;Xi=ai}\=j 



Note that the sum inside is a finite sum with at most (^) summands and it means that 
we sum over all those g{x) where there are exactly j arguments that come from the lower 
point a and the other ones come from b. As an example, one has for d = 1 the usual 
result g{A) = g{b) — g{a), for d = 2 one has 

9{A) =5(61,62) - g{hi,a2) -5(01,62) + 5(01,^2), 
and for (i = 3, ()2.ip becomes 

a{A) =g{hiMM) -5(^1,^2,03) -5(^1,^2,63) -5(^1,62,63) 
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+ g{bi,a2,as) +5(01,62,03) +5(01,02,^3) -5(01,02,03)- 

These formulas become much simpler if g is the product of d functions with one argument 
each, i.e. g{x) = YY^^i gi{xi). Then (|2.ip can be easily seen to reduce to 

d „ d 

9{^) = n / \a,,bj]{xj)g{dxj) = JJ {gj{bj) - gjiaj)), 
j=l j=l 

which is what one expects. This yields a measure fig on the Borel cj-field on by 
considering the Lebesgue extension of g in (|2.1|) . With respect to this measure one can 
now derive an integration theory for real valued functions in / G L^(//g) where p € [1, 00], 
and for such a function we denote the integral by 

/ f{x)g{dx) = / f{x)fig{dx), 
Jo Jo 

where O G Oc{U). This integration theory will be used in Sectional when we define a 
stochastic integral with respect to a random field, which has d+1 variables, d being the 
spatial dimension and one the temporal dimension. 

2.3 Hilbert-valued Malliavin Calculus 

In this section we provide some ideas and the main results we need related to Malliavin 
Calculus. It will however not be sufficient to only look at Malliavin Calculus for real- 
valued random variables or M*^- valued random vectors as treated extensively in [14] . 
Instead we have to deal with random elements taking values in some separable Hilbert 
space. Some sources for this are [H Section 5] or [TTj and references therein. A more 
general setting is the one treated in |16) where the authors treat Malliavin calculus for 
random elements taking values in some UMD Banach space. In this section we will 
without further notice identify the dual of a separable Hilbert space Q* with Q. 

Let ^ be a separable Hilbert space and let {W{h),h G Q) be an isonormal Gaussian 
process, see [HI Section 1.1.1] for some of its properties. Let be the probabil- 

ity space induced by the isonormal process. Furthermore we choose another separable 
Hilbert space Gi, and we consider the class of Hilbert-valued smooth random elements 
F G L'^in^Gi) given by F = f{W{hi),...,W{K)) for /ii, . . . , /i„ G ^, n G N and 
/ : M" — )• Qi which is infinitely Frechet differentiable with some boundedness condition, 
polynomially bounded or bounded. These functions / are dense in L'^{Q,;Qi). For these 
random elements F, the Malliavin derivative is given by 

DF := ^(W{hi), . . . , WiK)) hj. 
j=i ^1 

Consequently, for smooth random elements F we can interpret its Malliavin derivative as 
another random element with values in L2{Q, Gi), the space of Hilbert-Schmidt operators 
from G to Gi, or equivalently in the tensor product Gi^G since L2{G,Gi) is isomorphic 
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to this tensor product. Throughout this article we wih mainly work with the first 
approach, but occasionally use the second one when it is more convenient. We can 
also apply projections onto the coordinates of Qi leading to one-dimensional Malliavin 
calculus. In fact, for some / G Qi 

D'F := {DF, l)g, = (i^.iWihi), WiK)), l\ h,. 

A special choice for I in the previous equality is an element from a CONS of denoted 
by (efc)fceN yielding Dy^F := D'^'^F. A similar calculation for a CONS in Q leads to 
directional Malliavin derivatives as in the one-dimensional case. 

As in the real- valued Malliavin calculus, the operator D is closable in L^(i7;^i) and 
we define the Malliavin derivative of an element F G Lp'iyi; Qi) which can be represented 
as a limit of a sequence of smooth ^i-valued random elements (F„),k=n to be the limit 
of the Malliavin derivatives of the elements of the sequence, i.e. DF := lim„_>oo 
This convergence takes place in L'^{Q; L2{Q, Gi))- The space of all such elements will be 
denoted by D^'-^(Cyi) and it has the norm 

lli^lll2,e,=IE[||i^||^J+E[||Z)F||i^(^^^J, 

where the index Qi will be dropped if this does not cause any confusion. One can define, 
as in the real-valued case, the spaces B'^'P(^i) for /c G N and p > 1. 

We also need to define a Hilbert- valued equivalent to the divergence operator 6. This 
operator dg-^ : L'^{Q,; L2{Q,Gi)) — )■ L'^{Q;Qi) is defined to be the adjoint of D, i.e. 

E[{DF,G)L,^g^g^)]=E[{F,5g,{G))g,], 

for F G B^'^iGi) and ah G G L^{n; L2{G,gi)) for which 

\E[{DF,G)L,ig,g,)]\<G{B[\\F\\l])'/\ 

From now on we drop the index Qi from the divergence operator if this does not cause 
any confusion. 

Having defined these two operators we will now collect some calculus rules which we 
will rely on in the subsequent sections. First, we see that D and 6 are unbounded linear 
operators. This implies that one can pull bounded linear deterministic operators or 
functionals in and out of the Malliavin derivative and the divergence operator. As in 
the real-valued Malliavin calculus, there is a product and chain rule for the Malliavin 
derivative which in the Hilbert- valued case needs some explications. Let Qi, Q2 and 
be separable Hilbert spaces and let F G -^^2(^1)^2) and G G -^^2(^2, ^3) be two random 
linear operators which are Malliavin differentiable. Then GF G ][D^'^(L2((?i, ^3)) and 

D{GF) = {DG)F + GDF, (2.2) 

where this equality has to be interpreted as {D{GF))h = {DG)F{h) + G{DF){h) for ah 
h € Gi. A similar rule applies for directional Malliavin derivatives. For smooth Hilbert- 
valued random elements F and G this is shown in [11. Lemma 2.1] and the general 
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case follows by an approximation procedure by Hilbert-valued smooth random elements. 
The chain rule in Hilbert-valued Malliavin calculus is defined for functions (p : Gi ^ G2 
which are either Frechet differentiable or Lipschitz continuous. Let F € B^'^(^i). Then 
(t>{F) G ©^'2(^2) and 

D<p{F) = cp'{F)DF, (2.3) 

where (j) is Frechet diff'erentiable and (p' denotes the Frechet derivative of (p. If (j) is only 
Lipschitz continuous, then D(I){F) = ZDF where Z \s a random linear operator from Qi 
to Q2 whose norm is almost surely bounded by the Lipschitz constant of (j). 

Li our setting with a cylindrical Wiener process B on Qq, we can make some simplifi- 
cations. First we note that in this setting, the Hilbert space Q on which the isonormal 
Gaussian process is defined is equal to L2([0,T]; Qq) and one can reinterpret the Malli- 
avin derivative DF as some ^2(^/0) ^i)-valued stochastic process {F)tF)f^\f^ j'l^ on the time 
intervall [0, T] given by 

DtF := ^ ^{W{hi), WiK)) hj{t), 
j=i 

where hj € L'^{[0,T];Qo) for all 1 < j < n. Therefore DF actually denotes an equiv- 
alence class of functions from il. x [0,T] with values in L2{Gq,Gi), but one can find a 
representative such that DfF is measurable in lo for all t € [0, T] and that {DF){u}) is 
measurable in t for all cj E 17, which we denote as the Malliavin derivative of F. As in 
the general case, one can define the spaces n^'P{Qi), but moreover one can also define the 
spaces h'''P{gi) to be LP([0, T]; ©'=•^(^1)). In the classic real-valued Malliavin calculus 
L^'2 := L-'^'2(R). If = 1 then the norm is given by 

= rm\FtrgM'+ f f n\D,F,\\l^^^^^^^^]dsdt. (2.4) 

i/ t/ 

If > 2 then iterated Malliavin derivatives and further integrals are added to this 
expression. 

For the t?i -valued divergence operator 5g^ this has the consequence that it reduces 
to the t?i-valued Skorohod integral and for all G S L2([0, T] x Q.;L2{Qq.,Gi)) we write 
/o Gs5Bs instead of 5{G). If moreover G is predictable, then this integral turns out to 
be the ^i-valued Ito integral. 

The last issue we focus on in this section is the interplay between the Hilbert-valued 
Malliavin derivative and Skorohod integral. First, we have the general commutator 
relation D5{u) = u + 5{Du), similar to |14l (1.46)]. Let now u be a stochastic process 
in L-^'2(L2(t?0i ^1)) and we assume that for all t G [0,T] the process iDtu{s))s£[o^T] 
is Skorohod integrable and the process (J^ Dtu{s)6B{s))i;^[QT] has a version which is 
in L2(il X [0, T]; L2(^0) ^i))- This condition holds for instance if u is twice Malliavin 
differentiable. Then u{s)dB{s) G B'^'^iGi) and for ah t G [0,r] 

r-T rT 



Dt [ u{s)6B{s) = u{t) + [ Dtu{s)6B{s) 
Jo Jo 
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Finally, we provide a Hilbert- valued integration by parts formula which is inspired by 
[U Theorem 5.2]. However, we need it in a slightly more general setting which is why 
we include a quick proof here. Before we start, we fix a notation. Let throughout this 
article tr^p denote the trace of a linear operator A : Qq ^ -^2(^01 ^1) taken only over Qq, 
i.e. 

trgo(A) := ^{Aek,ek)go, 
km 

where {ek)km is a CONS of Gq. Consequently, the object trg(,(^) takes values in Qi. 
Prom the definition of the Skorohod integral as the adjoint of the Malliavin derivative 
one has that for all u G L'^{n x [0, T]; ^2(^0, ^1)) and A e B^'2(L2(ei, ^2)) 



E 



A [ u{s)6B{s) 
Jo 



E 



T 



trg, {{DsA)u{s))ds 



(2.5) 



where the integrand is ^2-valued and the integral is understood as a Bochner inte- 
gral. Similarly, one could write the trace outside the integral, which would yield an 
-^2(^0) Go <^ ^2)-valued integrand and integral. Now we are in the position to formulate 
the integration by parts formula. 

Proposition 2.2. Let A G D^'2(L2(^i, ^2)), and let ue L'^iQx [0, T]; L2(^o, ^1)) be in 
the domain of the Skorohod integral 5g-^^. Then Au € Dom((5g2) and 

[ Au(s)6B{s) = A [ u{s)6B{s) -trg^ [ Ds{A)u(s)ds, (2.6) 
Jo Jo Jo 

for all t G [0, T] . Note that under the conditions above the right-hand side of this equality 
is an element in L^(il;^2)- 



Proof. Assume that ^ G W'^{L2{Q2,'^))- We have to show that 



E 



^ [ Au{s)5B{s) =E "^(a [ u(s)5S(s) - trg^ [ Ds{A)u{s)ds 
Jo \ I \ Jo Jo 



(2.7) 



for all such ^. Then, a calculation similar to ^ Proposition 5.3] yields 



E 



^' [ Au{s)6B{s) 
Jo 



E 
E 
E 



tig, [ Ds{'^)Au{s)ds 
Jo 

iig, f [D,{-^A)--^D,A)u{s)ds 
Jo 

^{^^ u{s)5B{s) -txg^ D,Au{s)ds 



where we used 



in the first equality. This implies the assertion. 



□ 
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3 Stochastic integration 



In this section we provide an exact definition for the stochastic integral in ()1.2p with 
respect to an integrator as in (jl.ip . We are keen on deriving an integration theory that 
also covers singular g, i.e. where g{t, t) is not well-defined. In order to motivate the 
definition of the stochastic integral we provide a heuristic calculation that shows how 
each term comes into play. Throughout this and the following sections we work under 
the following assumption which have already been mentioned in Section [TJ 

Assumption 3.1. Fix T > and let t G [0,T], and let B = {Bt)t>o be a cylindrical 
Wiener process onl-Li. Furthermore, letg{t,s) G L{'H2,7i2) (nonrandom) for all <s< 
t <T and let {a{t))t>o be an L{T-Li,T-L2)-valued stochastic process such that g(t,s)a{s) £ 
L2{'Hi,T-L2) for all s G [0,t) and l^Qij(-)g{t,-)a(-) G Dom(5). Assume that for all s G 
[0,t) the L{'H2,'H2)-valued vector measure g{du,s) has bounded variation on [u,v] for all 
0<s<u<v<t. 

We note again that unlike in [3] we have not assumed that a is predictable with 
respect to -B, so that the integral in (jl.ip is a genuine Skorohod integral. In the following 
derivation we will first assume that 7^3 = M, the general case will be discussed shortly 
after. The basic idea for the calculations that follow is to expand all operators which 
appear in (jl.ip and (|1.2p into their coordinates, perform similar calculations as in [3] 
and then reassemble the original operators to get closed-form expressions for Hilbert- 
valued random elements. To this end we fix {ek)keN to be a CONS of and {fi)i£n 
to be a CONS of Tii. Then (i?');^^ := {{B, f 1)^,1) leN is a sequence of independent, 
one-dimensional Brownian motions. With the help of these two CONS, we can expand 
X'^ := {X, Ck) for all /c G N in the following way 



where in the last term {g(t, s),ek)'H2 denotes the linear functional from T-L2 to M which 
is defined by {g{t, s),ek)H2i^) •= {git, s)x,ek)H2 for all x G 'H2- The reason why we 
introduce this notation is to perform calculations in a more intuitive manner, in fact, 
this is a deterministic linear functional which commutes with the Skorohod integral and 
the Malliavin derivative as mentioned in Section [2.31 

Now we are in the position to motivate the definition of the stochastic integral. 
In what follows we first assume that the random integrand Y : [0,T] L(T-L2,T~{-3) 
is differentiable, but this assumption will be removed afterwards. In the following 
y(t)(efc) = Y{t)ek is just applying the linear operator Y{t) to G 7^2- In this case 
we obtain 




W2 lem 



= Y1 [\9it,s),ek)nMs){fi)SB\s), 



(3.1) 
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J2 (Y{t){ek)X\t) - I ^{u){ek)X\u)du 
fceN ^ JO ^ 

km k,im ^ " ^ 

j j Ds,i(^^iu){ek)^{g{u,s),ek)n2'^is){fi)dsdu 

Y,y{t){ek)x\t)- 

km k,im ^ ° 



-{u){ek){g{u,s),ek)H2(^{s){fi)SB^{s)du 



(3.2) 

where we first expanded the stochastic integral in ^2 along the coordinates of integrand 
and integrator, substituted and did a series expansion in Hi. Then we pulled the 
linear operator dY{u)/du inside the stochastic integral using Proposition [221 used the 
stochastic Fubini's Theorem, see |14l Exercise 3.2.7], pulled the deterministic bounded 
linear operator {g(t,s),ek) inside the Malliavin derivative and commuted the Malliavin 
derivative and the deterministic integral. Using similar steps and Proposition 12.21 again, 
we calculate the first term on the right-hand side of the last expression to be equal to 

Y,Y{t){ek)X''{t)= 5]y(t)(efc) f\g{t,s),ek)nMs)ifi)SB'{s) 
km k,im ° 

= E f fY{t){ek){g{t,s),ek)n2^{s){fi)5B'{s) 
k,im ^ 





Now we substitute this term into (|3.2p to obtain 
rt 



+ / D,^i{Y{t){ek)){g{t,s),ek)n:Ms){fi)d, 



[ Y{s)dX{s) 
Jo 

Q (Y{t){ek){g{t,s),ek)n2- ^(u)(efc)((7(tx, s), efc)^,du^ ct(s)(/0<55'(s) 



k,im 



* dY 

-Q:^iu)iek){giu, s),ek)-H2^u ) a{s){fi)ds 



+ Ds,i(Y{t){ek){g{t,s),ek)H2- 
E ( f (Y(^)iek){g{s,s),ek)n2+ /V(n)(efc)/|^(n,s),e,\ dn] a(s)(/0<5S'(s) 
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+ j Ds,i(Y{s){ek){g{s,s),ek)H2 + Yiu){ek)(^^iu,s),ek^ d'i?ja{s){fi)ds 
= 1^ (Y{s)g{s, s) + ^* Y{u)g{du, s^j a{s)6B{s) 

+ tvn, {Y{s)g{s, s) + ^* Y{u)g{du, s)^ a{s)ds, (3.3) 

by performing a deterministic integration by parts procedure using the fact that we can 
commute Frechet differentiation and the projection onto the k-ih. coordinate since they 
are bounded hnear operators which commute with the Frechet derivative. Then, we 
summed up over both CONS {ek)k&] and {fi)i^m. 

Next we briefly treat the case of a general separable Hilbert space ^{3. For this fix 
a CONS of ■7^3, denoted by {dk)k&%^ and use the expansion Y{t) = 'Ylke.'HO^i^)->'^k)dk- 
Then we perform all the calculations above on each coordinate of 7^3 separately and 
at the end sum up again to obtain a closed expression for the integral. Note that for 
this summing up to be true, certain summability conditions on the elements in each 
coordinate have to be assumed. 

These derivations motivate the definition of the following linear operator from 1^,2 to 
^3 

/Cg(/i)(t, s) := h{s)g{t, s)+ f {h{u) - h{s))g{du, s), (3.4) 

J s 

where the integral is defined as an integral with respect to an L('H2) ^2)-valued vector 
measure whenever it makes sense; note that g{du, s) has finite variation on all subin- 
tervals [v, t] where v > s hy definition. Sometimes we will call the operator in ()3.4p 
the kernel associated to X. Under the assumptions in 13.11 and the ones that follow in 
Definition 13.21 this linear operator is well-defined. We remark that this operator is the 
infinite dimensional analogon to the one which appears in [3], which in turn already 
appeared in [T]. In some special cases this operator can be rewritten. In fact, if g{s,s) 
is a well-defined linear operator from 7i2 to T-L2, then 

JCg{h){t,s) = h{s)g{s,s) + h{u)g{du,s), (3.5) 



Note that this is the kernel appearing in ()3.3p . and now it is obvious that ()3.4p is a 
generalization of it. If g{-, s) is absolutely continuous with respect to the one-dimensional 
Lebesgue measure on [0,t] with density (j){-,s), then 

]Cg{h){t, s) = h{s)g{t, s) + / [h{u) — h{s))(l){u, s)du, 

J s 

where the integral is understood as a Bochner integral. This situation applies in partic- 
ular if g{-, s) is Frechet differentiable. If 51 is a shift-kernel, i.e. g depends on t and s only 
through their difference, the kernel (|3.4p can be rewritten as 

/•t—s 

lCg{h){t,s) = h{s)g{t - s) + / h{u + s)g{du). 
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Going back to (j3.3p . we see that we can define the stochastic integral in (jl.2p as 

fY{s)dX{s):= f ICg{Y){t,s)a{s)6Bis)+tTn, f D,{]Cg{Y){t,s))a{s)ds. (3.6) 
Jo Jo Jo 

Next we describe the domain of this integral. 

Definition 3.2. Fix t > 0, let X he defined by (jl.ip and assume \3.1\ We say that 
a stochastic process (^(■s))sg[o,t] belongs to the domain of the stochastic integral with 
respect to X, if 

1. the process {Y{u) — Y[s))u^i^s,t] integrable with respect to g{du,s) almost surely, 

2. s /Cg(y)(t, s)(7(s)l[o,t] (s) is in the domain of the Hs-valued divergence operator 
5B, and 

3. JCg(Y){t, s) is Malliavin differentiahle with respect to Dg for all s € [0,t] and the 
T-L^-valued stochastic process s i-> tvy^^^ Ds[ICg{Y){t, s)^a{s) is Bochner integrable 
on [0, t] almost surely. 

We denote this by Y ^ l^{0,t) and the integral Jq y(s)dX(s) is defined by ()3.6p . 

Now one may ask about the the concrete form of the domain of this integral. So far, 
we have been unable to derive a characterization of it which is a similar problem as in 
anticipating calculus, where the domain of the divergence operator 6 is not completely 
characterized. Instead, one can identify a subset L^'^ in this domain. This is also the 
case for the X-integral that for some cases, we can identify a subset of its domain. In 
fact, if fj is assumed to be Malliavin differentiable then we can define the subset Xj'^2(0; 
which is given by the set of processes for which the seminorm 



\lCgiY){t,M-)h 



2,L2{Hi,H3) 



is finite. The equivalence classes of this seminorm depend on the exact shape of g and 
a. Similarly, one can define a semi-inner product in an obvious way. The setXi^2(0'*) 

IS 

then included in l^{0,t), since for some Y G Ij''2(0,t) we have 



E 



< 2E 



Y{s)dX{s) 



2 1 



Hz 



lCg{Y){t,s)a{s)5B{s] 



2 ■ 


+ 2rE 


[L 


«3- 







ixn, DslCg{Y){t,s)a{s)f^^ds 



<CT||/C,(y)(t,-)C7(-)||12;%' 



(3.7) 



where we have used the continuity of the Skorohod integral on L^'^('H3), see [HI (1.47)]. 
Note that the second term in the second line is a part of the Malliavin derivative of the 
integrand of the stochastic integral and therefore already included in the norm estimate 
for the Skorohod integral. This space Zj''2(0,t) takes the role of L^'^ from classical 
Malliavin calculus. 
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4 Calculus with respect to the integral 



In the first two subsections we present some general properties of the stochastic integral 
defined in the previous section. Afterwards, in subsections 14.31 and 14.41 we investigate 
some particular cases for the integrand. 

4.1 Basic calculus rules 

At first we can conclude that the integral defined in the previous section is linear. This 
follows immediately from the linearity of the Malliavin derivative, divergence opera- 
tor and Lebesgue(-Stieltjes) integrals. Formally we have for Y, Z ^ 1-^(0, t) and two 
constants a, 6 G M that aY + bZ e X^(0, t) and 

[ {aY{s) + bZ{s))dX{s) = a [ Y{s)dX{s) + b [ Z{s)dX{s). 
Jo Jo Jo 

Another immediate property follows from integrating constants. In fact, if we choose 
'H2 = ^3 and Y = id-^2) which is easily seen to be in l^{0,t) for all t > 0, then 
}Cg{Y){t, s) = g{t,s) for all s < t and since g is deterministic D}Cg(Y){t, s) = 0. This 
implies 

/ dX{s)= [ idn,dX{s)= f g{t,s)a{s)6B{s) = X{t). (4.1) 
Jo Jo Jo 

By combining this with the linearity of the integral for the deterministic integrand Y{s) = 

l[M,t)]id-H2 we obtain 



f l[„,„](s)dX(s) = / l[^^,]{s)idn,dX{s) = X{v)-X{u) 
Jo 



(4.2) 



which gives us the intuitive property that the integral over an indicator function is the 
increment of the integrator. We have furthermore that if < t < T and Y G 1-^(0, t) 
then yi[o,4] GX^(0,r) and 

r l[o,t]is)Yis)dXis)= /V(s)dX(s). (4.3) 
Jo Jo 

This can be seen by splitting the integral jj--- g{du, s) into ■ ■ ■ g{du, s) and ■ ■ ■ g{du, s) 
Note that this equality does not hold if we only assume that Y G 1-^(0, T) because the 
fact that X is Skorohod integrable over [0, T] does not in general imply that -'^Ifct] is 
Skorohod integrable over [0,T], see [14^ Exercise 3.2.1]. However, if a is assumed to be 
Malliavin differentiable, this does hold. In fact, in this case JCg{Y){T, s)a{s) is Malliavin 
differentiable for all s G [0,t] which implies that JCg{Y){t, s)a{s) is Skorohod integrable 
over [0,t] for all t G [0, T]. Combining (j4.3p with the linearity of the integral we have 
immediately that for < u < v < t and Y G Z^(0,n) nl^{0,v), that yi[n,t,] G l^{0,t) 
and ^ 

/ y(5)l[,,,](s)dX(s)= rYis)dxis)- ry(s)dx(s). 
Jo Jo Jo 
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Using these basic rules we can derive more interesting properties for the stochastic 
integral with respect to X. 

Proposition 4.1. Assume that Assumption \3.1\ holds and that X is defined by (jl.ip . 
Let t > and assume Y G l'^{0,t). 

1. Let Z he a random linear operator from T-L^ to another separable Hilhert space T-L^ 
which is almost surely hounded. Then ZY G X^(0, t) and 

[ ZY{s)dX{s) = Z [ Y{s)dX{s) almost surely. 
Jo Jo 

2. The X -integral is local, i.e. ifY{s) = for allO < s <t almost surely, then 

/ Y{s)dX{s) = almost surely. 
Jo 



3. Let Y he a simple process, i.e. Y = ^jl(tj,tj+i] where Zj is a random linear 

operator from 7^2 to which is almost surely hounded for all j = 1, . . . ,n — l and 
< ti < . . . < tn ^ t is a partition of the interval [0, t]. Then Y G T^{0, t) and 



„t n-l 

/ Y{s)dX{s) = Y,Zj{X{t,+,)-X{t,)). 



4. Let furthermore a he Malliavin differentiable. Then the X -integral is a continuous 
linear operator from X^2{^^'t) to L'^{Q.]'H^). 

Proof, (i): Note that by definition of Kg we have K,g{ZY){t, s) = ZK.g{Y){t,s) and by 
the product rule of Malliavin calculus ()2.2|) 



Ds{Z}Cg{Y){t, s)) = Ds{Z)]Cg{Y){t, s) + ZDslCg{Y){t, s). 
This and the Hilbert- valued integration by parts formula ()2.6p yield the assertion 

[ ZY{s)dX{s) 
Jo 

= I Z)Cg{Y){t,s)a{s)SB{s)+tvn, [ {ZDs}Cg{Y){t, s) + Ds{Z)}Cg{Y){t, s))a{s)ds 
Jo Jo 

= zQ ICg{Y)it,s)a{s)6B{s)+tTn, 1^ DsJCg{Y){t, s)ais)ds^ 

= Z I Y{s)dX{s). 
Jo 

(ii): This claim follows immediately from the fact that the Malliavin derivative, the 
Skorohod integral and the Lebesgue(-Stieltjes) integral are local operators. 

(iii) : This is a combination of (i) and (j4.2|) . 

(iv) : This is an obvious conclusion from (i) and the continuity in (|3.7p . □ 
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Another property we want to investigate in this section are two projection equahties 
with respect to the CONS {ek)k€N of 'H2- The first one is given by the apphcation 
ICg{Y)(t, s) to Cfc, which yields 

)C,iY){t,s)iek) = ICg^edY){t,s) (4.4) 

which is easy to see from (|3.4p . where g{ek) G ^2 is defined by g{t, s)ek for all < s < t < 
T. The other projection equality is integration with respect to = {X, efc)-^2) which is 
a real-valued stochastic process. Applying the definition of the stochastic integral (|3.6p 
invoking ()3.ip with = 1^ we see 

\{s)dX\s)= f )C^g^,^){Y){t,s)SB{s)+tvn, f DsK.^g,e,){y){t,s)ds, 
Jo Jo 

where y is a stochastic process taking values in a Hilbert space H^. Note that if Ti-^ = 
M, the function s 1— )• /C^^ g^) (^)(^) ■s) is real-valued, so we can relate this to the one- 
dimensional case treated in ^J. In order to recover exactly that situation, one has to 
assume Tii = M, too. 



4.2 Semimartingale condition 

In this section we focus on stating a condition under which the stochastic integral process 
1 1—)- X{t) is a semimartingale. This condition is a regularity and smoothness assumption 
on g. The following proposition is inspired by [21 Proposition 5]. 

Proposition 4.2. Let t > and assume that g{t,s) is well-defined for all < s < t. 

Furthermore, assume that there is a hi-measurahle function (p : [0,T] — t- L{7i2,'H2) such 
that 

9it,s) = gis,s) + (p{v,s)dv, (4.5) 



for all < s <t, where this integral is defined in the sense of Bochner and 

rt ru 

\9is,s)\\lrn,^^y^.ds<oo and / / {{(piu, s)\\lrn,^^^.dsdu < 00. 



Suppose furthermore that a is adapted to B and pathwise locally hounded almost surely. 
Then X defined by (|l.ip is a semimartingale with decomposition 



X{t) = f g{s,s)a{s)dB{s)+ f f 
Jo Jo Jo 



t rs 

cl}{s,u)a{u)dB{u)ds. (4.6) 

Jo 

Proof. Substituting the definition of g and applying a stochastic Fubini Theorem we 
obtain 

/*i /*i /* s 

X{t)= / g{t,s)a{s)dB{s) = / g{s, s)a{s)dB{s) + / / (l){s,u)a{s)dB{u)ds, 
Jo Jo Jo Jo 

which can be seen to be a semimartingale by \13\ Remark 26.4]. □ 
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Prom ()4.5p we see that (j) can be interpreted as some sort of derivative oi g. In fact, if 
g{-,s) is Frechet differentiable, then 4>{t,s) = ^{t,s). Note that one can also relax the 
condition on a to be adapted, leading to so-called Skorohod semimartingales, see [15] . 

In the case when X is a semimartingale, the question arises in which relationship our 
integral and the classical integral with respect to the semimartingale X, see [17j, stand 
to each other. We suppose in the following that Y is predictable, which is a requirement 
in the integration theory with respect to semimartingales. Using the decomposition of 
X given in (|4.6p . then the stochastic integral of Y with respect to X, denoted hy Y ■ X, 
is given by 



{Y-X)t= / Y{s)g{s,s)ais)dB{s)+ / Y{s) / (p{s,u)a{u)dB{u)ds 
Jo Jo Jo 

However, by applying integration by parts and the stochastic Fubini Theorem, we get, 

t f-t f-t 

Y{s)gis,s)a{s)dB{s)+ / / Y{u)^{u, s)dua{s)6B{s) 

Jo Js 

+ tr^j / / DsiY{u))(l){u,s)a{s)duds 

Jo Js 

= [ Y{s)dX{s)-trn, [ Ds{Y{s))g{s,s)a{s)ds. (4.7) 
Jo Jo 

Therefore, these integrals are equal only if DsY{s) = for almost all s € [0, t]. This holds 
for instance when y is a simple process as in Proposition [530 where Zt- is assumed to be 
measurable with respect to iBt)t£[o,ti]- Similarly, this also holds when Y is deterministic, 
or independent of the driving Wiener process or lags behind the driving Wiener process. 

The difference of the integral comes from the following fact. Starting from the pre- 
dictable simple processes for which DsY{s) = for almost all s G [0,t], we define the 
semimartingale integral as the closure of the simple processes under the L^([X]'^)-norm. 
For the X-integral, this is not what one does. In fact, predictable simple functions do not 
even approximate all the processes in the domain of the X-integral with respect to the 
L^([X]'^)-norm. For instance, take Hi = 'H2 = "Hs = M, let y be a Brownian motion on 
[0,T] and let 5 = 1 and a = 1. In this setup, ICg(Y){t, s) = B{s) and DslCg{Y){t, s) = 1. 
Therefore, Y G 1^(0, t) for ah t G [0,r] and 



/"V(s)dX(s)= fB{s)dB{s)+ f ds = l{B{tf +t). 
Jo Jo Jo ^ 



(4.8) 



But the Malliavin derivative of each predictable simple process Yn in Proposition I4.1IBI 
approximating the Brownian motion is equal to zero and therefore, does not converge 
to 1. Another set of integrands on which the integrals conincide are semimartingales of 
the form 



Y{t) = At+ [ u{s)dB{s), 
Jo 



where A is a process of bounded variation and we assume that DsA(s) = and Dsu{s) 
for almost all s G [0,T]. 
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4.3 Deterministic Integrands and OU processes 

In this section we investigate the integral Jq h{t, s)dX{s) if h is chosen to be deterministic. 
To this end fix t > and let s i— >• h{t, s) be a measurable function from [0, t] into 
L{H2-,T~i^) such that h{t,u) — h{t,s) is integrable with respect to g{du,s) over [s,t] and 
Definition [230 holds. In this case the kernel K,g{h) is deterministic and has therefore a 
Malliavin derivative equal to zero, which immediately implies [3] in Definition 13.21 Then, 
by the definition of the integral we have 

t f-t 
h{t,s)dX{s)= / ICg{h{t,-)){t,s)a{s)6B{s) 
Jo 

h{t,s)g{t,s) + {h{t,u) - h{t,s))g{du,s)^a{s)6B{s). 

Note that this is a new process of the same type as X in (jl.ip where now ICg{h) plays 
the role of g. 

Now we apply this to a specific deterministic integrand which is stationary in time in 
order to show that Ornstein-Uhlenbeck (OU) processes driven by X are representable as 
a stochastic integral of a deterministic integrand with respect to X. By an OU process 
we mean the solution {Y{t))t>o to the infinite dimensional SDE dY(t) = —AY{t)dt + 
FdX(t), or in integral terms 

Y{t) = - [ AY{s)ds + FX{t), (4.9) 
Jo 

where {X{t))t>o is the process defined in ()l.ip and —A is an unbounded linear operator 
from 7^3 to whose domain is dense in Ti^ and that generates a Co-semigroup on 
denoted by S{t) = e~^^ for all t > 0. In the previous equation F is a bounded linear 
operator from to H^. Note that this is a first step into the investigation of SDE 
driven by a Volterra-type process. As suggested from the classic situation where the 
driving noise in (j4.9p is a Wiener process we obtain that the solution to (j4.9p is given 
by a stochastic convolution. 

Proposition 4.3. Suppose that Assumption \3.1\ holds. Assume furthermore that t i— 
a{t) is pathwise locally bounded almost surely. Assume for all s G [0, T] that the map 
(s,T) 9 u I— )• (e"*^""*)"^ — e"^"^) is integrable with respect to g{du,s), that 

2 

FJCg{S){s,u)a{u)5B{u) 



E 

and 



ds 

Hz . 



< oo, (4.10) 



lim II {idu, - e-(*~^)^) <7(t, s)\\^^ = 0. (4.11) 



Then S G l^{0,t) and 

Y{t)= f S{t-s)FdX{s)= fe-^'-'^'^FdXis) (4.12) 
Jo Jo 



solves the SDE (US 
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First let us explain the two conditions ()4.10p and ()4.1ip . The former one is the Hilbert 
valued equivalent of the one in [14. Exercise 3.2.7] which is sufficient for an application 
of the stochastic Fubini theorem with the divergence operator that we will use in the 
following proof. If a is assumed to be predictable, the Skorohod integral becomes an Ito 
one since ICg{h) is deterministic and in this case one can derive similar as in [3l, Lemma 
2] that 

T rT 

^ duds < oo 



is a sufficient condition for the stochastic Fubini theorem to hold. The latter condition 
()4.1ip can be interpreted as a kind of weighted Co-semigroup condition. In fact, we 
demand that the semigroup S{t) = e~^^ goes to the identity as t i in a family of 
norms with weights given by the functions g{-,s). 

Proof of Proposition \4-^ Let throughout the proof G Dom(^*). By substituting (|4.12p 
and the definition of the stochastic integral into (|4.9p and applying the stochastic Fubini 
theorem one obtains 

C, / AY{u)dv) = I (a*C, [ }Cg{SF){u,s)a{s)6B{s)) du 
Jo / Jo \ Jo I n-) 



Hz JO \ JO I H2 

ft \ * 

AlCg{SF){u,s)du) C,a{s)5B{s) 



Hi 



and similarly for the first term of 



{C,Y{t))n, = l\}Cg{SF){t,srC,a{s)6B{s))n, 



These two equalities yield that 



C, Y{t)+ f AY{u)dvS = f I (lCg{SF){t, s)+ f AlCg{SF){u, s)dv\ C, (j{s)5B{s) 

Jo / Us Jo \ \ Js J / H2 

(4.13) 

Now we investigate the term in the brackets in the previous equality. Written out, this 
term is equal to 



e-(*-«MF(7(t, s) + j |^e~(*-")^ - e'^*"^)^ j Fg{du, s) 

+ j^^ (^Ae-("-")^F5(n, s) + (^e-("-^)^ - Afi-'^'^-'^^)Fg{dv, s)^ du. (4.14) 
Using Fubini's Theorem and [10', Lemma II.1.3(ii)] we see that 
f r(^e-("-^)^ - Ae-("-")^)F5(dv,s)dn 

J S J S 

- ^e-(^-^)^ + ^e-("-^)^ duFg{dv, s) 
au du I 
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= - f (e-(*-")^ - e-^'-''^^)Fg{dv, s) + [\l - e-^''-'^'^)Fg{dv, s). (4.15) 

J s J s 

Note that by an integration by parts procedure and using [101 Lemma II.1.3(ii)] again, 
the second term on the right-hand side of the last equahty may be calculated to be 

r (1 - e-(^-^)^)F(7(d7;, s) = {l - e-^'-''>^)Fg{t, s) - f Ae~^''-'^^Fg{v, s)ds. (4.16) 

J s J s 

where we used the condition ()4.1ip so that the boundary term at s does not appear. 
Assembling ()4.14p - ()4.16p we obtain 

}CgiSF){t,s) + f AJCg{SF){s,u)ds = Fg{t,s), 

J u 

which in turn, by substituting this into ()4.13p and using the fact that F is bounded, 
implies the assertion 

C, Y(t) + r AY(u)du) = /c, r Fg{t, s)a{s)6Bis)) = {(, FX{t))ns ■ 
Jo I \ Jo I 



The proof ends by the fact that A and hence A* are densely defined in which implies 

dm □ 

Note that (j4.12p has the exact same form as in the classic case when the driving process 
is a Wiener process. We however have to remark that this does not show uniqueness of 
the solution to the SDE (lOD. 



4.4 Volterra Processes as integrands 

In this section we want to calculate the stochastic integral with respect to X when 
the integrand is a Volterra process itself. This will also lead us to an integral of type 

X{s)dX{s). Throughout this section we assume that a = id-^j ^'^d assume that 
S is a (cylindrical) Wiener process on 7^2- However, one can reintroduce a into the 
calculations by formally considering adB as the noise. Moreover, we assume in the 
following without loss of generality that = M. The case of general separable Hilbert 
spaces then follows by the expansion Y{s) = ^^(^(s), /fc)/fc, where {fk)keN is a CONS 
of Ti^, and performing the subsequent calculations coordinatewise. 

Let in the following h be an L(?^2) -^(^2) IR))-valued deterministic function in t) 
and by the results of the previous section 

Y{s)= [ h{u)dX{u)= [ lCg{h){s,v)dB{v). (4.17) 
Jo Jo 

Here the integral is in fact a classic Ito integral. Furthermore, let {ek)keN be a CONS of 
T-L2- Using the definition of the stochastic integral with respect to X and the notion of 
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the projection in ()4.4p . we have 

/"V(s)dX(s) = V f f K:g{Y){t,s){em)5B^{s)+ [' Ds,„,JCg{Y){t,s){em)ds]. 

(4.18) 

Now, taking advantage of the specific form of Y, we calculate the integrand in the 
stochastic integral in (|4.18|) to be 

= Y{s)g{t,s){en,)+ [ {Y{u) - Y{s))g{du, s){em) 

J s 

}Cy{h){t,s)dB{v)g{t,s){e^)+ I I {)Cg{h){u,v) - }Cg{h){s,v))dB{v)g{du, s){em] 

Js Jo 

+ ICg{h){u,v)dB{v)g{du,s){em) 

J S J S 

I }Cg^eoW{t,s)g{t,s){e^)dB'{v)+ f I K:gi^,^){h){u,v)g{du, s){e^)dB\v) 



+ / / {^g{ei){h){u,v) - Kg(ei){h){s,v))g{du,s){eTn)dB {v)]. 



Js 



Y.{ f ^a{e^)^^9(e,)m;v)){t,s)dB\v)+ [ [ Kgi^,^^{h){u,v)g{du, s){e^)dB\v) 
E f }C^g'^ms,v,t)dB^(v), (4.19) 



where we have set 

^grg2(^)('5'^'*) = ^{v<s}K^g2iem)(.^gi{ei)ih){-,v)){t,s) + l{^y,} j K.g^i^^^){h){u,v)g2{du, s){em) ■ 

J V 

(4.20) 

Taking the Malliavin derivative in (j4.19p yields 



D,,^ICg^,^^{Y){t,s) = DsJYI f IC'g^{h){s,v,t)dB' 



(v)] =lCl^^{h){s,s,t). (4.21) 



Note that on the right-hand side the argument (s, s, t) appears which implies the explicit 
appearance of g{s, s) in the formulas, see below for a comment. Assembling all the terms 
- ()i:2T]l we obtain 

t 

Y{s)dX{s) 

E ( fH f iC'g':^{h){s,v,t)dB\v)5B^{s)+ f ^^{h){s,s,t)ds 
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= T. I I i<:''r^h){s,v,t)AB\v)6B^{s) + Y, [ IC^^ms, s,t)ds 

= E f [\ic's':^ms,v,t)+}c^j{h){v,s,t))dB^{v)dB^{s) + Y. fKrms.s^m 

If we now define the linear operator hy iCg^g{h){s,v,t){ei(S>em) '■= ^g^{h){s,v,t), 

we see in total that 

t 

Y{s)dX{s) 

t fS rt 

{iCg,g{h)(s, V, t) + iCg,g{h){v, s, t)) dB (v) dB (s) + tr^, / }Cg,gih){s, s, t)ds. 

Jo 

(4.22) 







JO 



The most important observation from this calculation is that every element in this 
formula is deterministic and therefore we see that the integral Jq Y{s)dX{s) belongs to 
the second Wiener chaos plus a term from the zeroth Wiener chaos. Moreover we see 
that the term )Cg^g{h)(s, s, t) appears in the correction term of the formula. This implies 
that the integral term in ()4.20p does not appear in this term, but in the first summand in 
(j4.20p the terms g{s, s) will appear explicitly, which is why we have to assume sufficient 
regularity for the kernel g. The reason for X to be regular here is that it appears in the 
integrand ()4.17p as well as in the integrator. Otherwise both would in general have a 
singularity at the same point which would cause problems when integrating. 

Similarly to these calculations one can easily extend the results to iterated X integrals 
of higher order. As long as a is deterministic, the X-integral of an integrand Y , which is 
an iterated X-integral, consists of a Wiener integral over lCg(Y) and a pathwise integral 
of the Malliavin derivative of ICgiY). This increases the number of the Wiener chaos of 
the integrand by one (Wiener integral) and decreases the number of the Wiener chaos of 
the integrand by one (pathwise integral). Therefore, for k odd (even), the A;-th iterated 
X-integral of a deterministic function is an element of all odd (even and zeroth) Wiener 
chaoses up to order /c, respectively. 

Now we apply this to the problem of calculating X{s)dX{s). Since both elements 
are 7^2-valued, we cannot apply ()4.22p verbatim. However, by Riesz Representation 
Theorem, we interpret the integrand as a linear functional on 7^2 and therefore 

f X{s)dX{s) := f X{srdX{s)= f\x{s),dX{s)) = Y, f X\s)dX\s). 

JO JO Jo A;GN-^° 

In this case we have h{s) = id-^j which implies ICg{h)(t, s) = g{t,s). We assume here 
that g is regular. Since X^{s) = j^{g{s,u),ek)'^,.^dB{s), we obtain that Dg^iX^^u) = 
{g{u,s)ei,ek)H2- Then, 
t-t 

^ ' X^{s)dX''{s) 



fcGN 



'0 
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+ E f Ds,i{x\s){g{s,s),ek)H.+ f X''{r){g{dr,s),ek)n?j(.ei)ds 
[ ( [ {g{s,u)dB{u),g{s,s)dB{s))n,+ [ [ ([ g(,r,u)dB{u), g{dr, s)dB{s: 

Jo \ Jo Jo Js \Jo I Ui 

+ ^^H2 { {g{s,s),g{s,s))u2 + (5(^s),5'(dr,s))^2 Ids. (4.23) 



We want to find a formula that links this term with X(t))-^2 ) which can be cal- 

culated to be 

^{X{t),X{t))n2 

= 11 {git,u)dB{u),g{t,s)dB{s))n2 
Jo Jo 

= [ [ {g{s,u)dB{u),g{s,s)dB{s))n2+ [ III g{r,u)dB{u),g{dr,s)dB{s)) 

Jo Jo Jo Js \Jo I Ui 

+ f r ( f g{dv,u)dB{u),g{s,s)dB{s)) - f f ( F g{t,u)dB{u),g{dr,s)dB{s) 

Jo Jo \ Js I ■H2 -^0 Js \ Js 

+ t f ( [\9{t,u) - g{r,u))dB{u),g{dr,s)dB{s)) (4.24) 

Jo Js \Jo I Hi 

using the order u < s < r < t. Note that the first and second term on the right-hand side 
of this equality are equal to first and second term on the right-hand side of ()4.23p . The 
three last terms in (|4.24|) can be rewritten by grouping the first and third one together 
as 




/ g{dv,u)dB{u),g{s,s)dB{s)) + [ ( [ g{dv,u)dB{u), g{dr, s)dB{s] 

10 Jo \ \ J s I Ui Js \Jr I Ui 

- f f I r g{r,u)dB{u),g{dr,s)dB{s) 

Jo Js \ J s I T-ii 

which we can calculate to be, using Fubini's theorem 

ft 



X] \\ g{dv,u)ei,ek) {g{s, s)em,ek)H2 

fc,/,mGN-^0 Jo \\Js /h2 

^ Is i.lr (5'(dr,s)e„,efe)^2^d5'(n)d5'"(s) 
{g{r, u)ei,ek)H2dB\u){g{dr, s)em,ek)H2dB'^is) 



k,LmeN ' 
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j[ [ '^ioieiMn, ( / iaidv, n)em, ek)n,^ {t, s)dB\u)dB^ (s) 



~ 12 {g{r,u)ei,ek)n2{9{dr,s)em,ek)H2dB^{u)dB"'{s) 
- f f [\g{r,u),g{dr,s))n2dB{u)(S)dB{s), 

Jo J s Ju 

where g{-, s)* means the adjoint hnear operator of g{-, s). The third term in ()4.23p which 
can be seen to be equal to 

XI / [{9is,s)ei,ek)H2{9{s,s)ei,ek)n2+ (fi'l^ s)ei, efc)^2 (^(dr, s)e/, efc)^2 j ds 
Yl [ ^{g{ei),ek)n2i(9i-:S)ei,ek)'H2){t,s)ds (4.25) 



k,l€N 

trw2 / ^g(5(-,s)*)(t,s)ds 



One could also reformulate the the integral term in the first line of ()4.25p . which has the 
form f{u)f{du) = - f{s)), to be 

{ {9{s,s)ei,ek)H2{9is,s)ei,ek)H2+ / ■s)e/, efc)^^ (^(dr, s)ei, 6/^)^2 I ds 

1 /■* 

= 2/^ i\\9it,s)\\l^(m) + \\9is,s)\\l^im))^'- 
So, collecting all the terms above we finally have an equality 

1^ X{s)dX{s) = ^ {X{t),X{t))n2 - 1^ ]Cg I' g{dv, u)*^ {t, s)dB{u) ® dB{s) 

- f f [\9{r,u),gidr,s))n2dB{u)0dB{s) + trH2 T /Cg((7(-, s)*)(t, s)ds, 

Jo Js Ju Jo 

which gives us the desired link between the X-integral of X and 

5 An ltd Formula 

In this section we derive Ito formulas for the processes X and Z defined in (jl.ip and 
()1.2p . In order to do this we rely on the Ito formula in Hilbert spaces with anticipating 
integrands in [11', Proposition 4.10]. 
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Proposition 5.1. Let T-L,JC be separable Hilhert spaces and B be a T-i-valued cylindrical 
Wiener process. Let moreover F € C^(/C;M) (the twice Frechet differentiable Junctionals) 
and let (V'(t))t>o be the stochastic process defined by 



V{t) = V{0)+ [ A{s)ds+ [ C{s)SB{s), 
Jo Jo 



lo Jo 

where V{0) G B^'^{JC), A G L1'4(/C) and C G L2'P(?^,/C) for some p > A, see Then 



F{V{t)) = F{V{0))+ f F'{V{s))A{s)ds+ f F'{V{s))C{s)6B{s) 
Jo Jo 

+ i trn 1^ F"iVmD-V)is){Cis))ds 



+ -tvn F"{V{s)){C{s)){C{s))ds 

where 

{D-V){s) = 2DsV{0) + 2 1 DsA{r)dr + 2 1 D,C{r)6B{r) 

Jo Jo 

Note the definition of the terms in the above formula. Since F' G -L(/C, L(/C, M)) the 
first two integral terms take values in M, as F(y{0)) does. For the forth term on the right- 
hand side we use the definition of F" to be an element of L{IC, L(IC, L{IC, M))) whereas for 
the fifth term we use the equivalent formulation F" G L{IC ^ /C, L2(/C, M)). Then, since 
{D-V){s) G L2{n,lC) and C{s) G L2{H,iq, we have that F"{V{s)){D-V){s){C{s)) G 
L2{n,L2{'H,R)) and F"{V{s)){D-V){s){C{s)) G L2{n, L2{'H,R)) so that the trace over 
Ti is in both cases well defined. Note that this trace is our way of writing the scalar 
products in the original formula in |1H Proposition 4.10]. 

The first issue in this section is to extend the above proposition to functions F G 
C^(/C; ICi) where /Ci is another separable Hilbert space, not necessarily M. We can reduce 
this to applving Proposition lS.ll coordinatewise in the following sense: set F'^ := (F, efc)/Ci 
and observe 

FHV{t)) 

= F'{V{0))+ [\F'y{V{s))A{s)ds+ [\F'yiV{s))C{s)dBis) 



+ -tvn {F^nV{s)){D-V){s)C{s)ds + -i.n {F'^)"{V{s)){C{s))C{s)ds 
= /f(F(0))+ f F'{V{s))A{s)ds+ f F'{V{s))C{s)5B{s) 



+ \ trn F"{Vis)){D-V)is)Cis)ds + ^ tr« F"iVis))iC{s))C{s)ds, e^ 
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This holds true since one can commute Frechet derivatives and projections, as they are 
bounded hnear operators. Hence we have identified each coordinate and summing over 
these coordinates yields an Ito formula in the Hilbert-valued case. 

Now we derive an Ito formula for processes like Z in ()1.2p which may be Hilbert-valued 
or real- valued. In order to apply Proposition 15. 1 1 to our case we will have to assume that 
g{s,s) is a well-defined linear operator for all s € [0,T] and that there exists a function 
(p as in (|4.5p . Therefore, the Ito formula in this section will only hold for the case when 
X is a Skorohod semimartingale, that is, with g satisfying the conditions of Proposition 
4.21 but with an anticipating integrand in the stochastic integral. So we will apply the 
above proposition to the following stochastic process 

t 

Y{s)dX{s) 

]Cg{Y){t,s)a{s)5B{s)+ f tin, {D,lCg{Y){t,s)a{s))ds 

Jo 

Y{s)g{s,s)a{s)6B{s) + [ [ Y{u)^{u, s)dua{s)6B{s) 

Jo Js 

+ / trn,D,{Y{s))g{s,s)a{s)ds+ / tr^, / DuiY{u))<j){u, s)dua{s)ds 

Jo Jo Js 

I b{s)5B{s)+ / ( / a(s,n)(5i?(u) +tr-Hi (6d(s)) +tr^i / aD{s,u)du]ds, 
Jo \Jo Jo J 







where we have used the stochastic Fubini Theorem and 

6(s) = Y{s)g{s, s)cr(s), a{s, u) = Y(s)(l){s, u)a{u), 

6z)(s) = Ds{Y{s))g{s, s)a{s), aois, u) = Ds(Y{s))(j){s, u)a{u). 

So we apply 15.11 to the following processes 



C{s) 



b{s) and A{s) = / a{s, u)SB{u) + ti-^-^ b£){s) + ti-^-^ / a£){s,u)du. (5.1) 
Jo Jo 



After the following theorem which sets up an Ito formula for the integral with respect to 
X, we will provide some sufficient conditions so that these processes satisfy the conditions 
in Proposition 15.11 



Theorem 5.2. Let T-i^ he a separable Hilbert space and let F : Tis 71^ he twice 
Frechet differentiahle. Let X and Z he defined as in (jl.ip and ()1.2p where we suppose 
that Assumption \3.1\ holds. Furthermore assume that g satisfies the semimartingale 
conditions in proposition \4-S\ Assume that Y and a are twice Malliavin differentiahle 
and 

C{s) eh^'PinuHs) and A{s) Gh^'^ns) (5.2) 
where A{s) and C{s) are as in dSTj. Then F'{Z)Y G l^{0,t) for all t G [0,r] and 

F{Z{t)) = F(0) + f lCg{F'{Z)Y){t,s)a{s)5B{s) +tvH^ f Ds{JCg{F'{Z)Y){t,s))a{s)ds 
Jo Jo 
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- ^ tru, 1^ F"{Z{s)) {Y{s)g{s, s)a{s)) {Y{s)g{s, s)a{s))ds (5.3) 

F(0)+ f F'{Z{s))Y{s)dX{s) 
Jo 



^tvn, I F"{Z{s)){Y{s)g{s,s)a{s)){Y{s)g{s,s)a{s))ds 



ft 

2 ^ 

Proof. Applying Proposition 15.11 yields 



F{Z{t)) =F{0)+ / F'{Z{s))b{s)6B{s) + / F'{Z{s)) / a(s, u)5S(u)ds 
Jo Jo Jo 

/*t ft/ f S 

+ / F'{Z{s))ixn,hD{s)ds+ / F'{Z{s))iTn, / az3(s,n)dnds 
JO JO Jo 

+ iiHi I F"{Z{s)) I Dsb{u)6B{u)b{s)ds 

Jo Jo 

+ tr^, / F"{Z{s)) I Ds I a{u,r)5B{r)dub{s)ds (5.4) 

Jo Jo Jo 

+ ti:m I F"{Z{s)) I DsiTH,hD{u)dub{s)ds 
Jo Jo 

pt ps pu 

+ / F"{Z{s)) / D.trn, / aDiu,r)drdub{s)ds 
Jo Jo Jo 

+ ^tTn, j'^F"{Z{s))b{s)b{s)ds. 
Note that the third term on the right-hand side of (|5.4|) can be rewritten using (|2.6|) as 

F'{Z{s)) I a{s,u)5B{u)ds 
Jo 

F'{Z{s))a{s,u)5B{u)ds+ [ [ trui {{DuF' {Z{s)))a{s,u))duds. 

Jo Jo 

(5.5) 



t ps 

Jo 



By the Ito formula in Proposition 15. II we know that all the terms on the right-hand side 
of ()5.4p are well-defined. Now we calculate the terms on the right-hand side of ()5.3p 
assuming that they exist. We will see that they are equal to a sum of terms on the 
right-hand side of ()5.4p . and we already know that the latter terms are well-defined. 
Therefore we can conclude that all the terms in ()5.3p are also well-defined and that the 
equality in ()5.3p holds. In fact, we have by definition 



/*t ft ft fS 

I }Cg{F'{Z)Y){t,s)a{s)dB{s) = / F' {Z{s))b{s)5B{s)+ / F' {Z{s))a{s,u)5B{u)ds 
Jo Jo Jo 

(5.6) 



and 



tin, f {DslCg{F'{Z)Y){t,s))a{s)ds 
Jo 
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tin, (^DsF'{Z{s))b{s)ds + F'{Z{s))bD{s) + I' {{DsF'{Z{u)))a{u, s) 



+ F'{Z{u))aDiu, s))dujds 



tin, [ {DsF'{Z{s)))b{s)ds + trn, I F'{Z{s))bD{s)ds 
Jo Jo 

ft f S ftfS 

+ trH, / / {DuF'{Z{s)))a{s,u)duds + tiu, / / F'{Z{s))aDis,u)duds. (5.7) 
Jo Jo Jo Jo 

Note that by the chain rule (|2.3p 



b{s)+ / D,b{u)6B{u)+trni / / a{u,r)6B{r)du 
Jo Jo Jo 

ps ps pu \ 

+ / DsticnibD{u)du+ / Dsti-Hi / aD{u,r)drdu]. 
Jo Jo Jo J 

With this we can rewrite the first term on the right-hand side of (j5.7p to 

tr«, C [D,F\Z{s)))b{s)ds 
Jo 

ft ft f S 

= tTHi / F"{Z{s))b{s)b{s)ds + tTm / F"{Z{s)) / Dsb{u)6B{u)b{s)ds 
Jo Jo Jo 

pt ps pu 

+ trni / F"{Z{s)) Ds a{u,r)5B{r)dub{s)ds (5.8) 
Jo Jo Jo 

+ trm [ F"{Z{s)) [ D,tTn,bD{u)dub{s)ds 
Jo Jo 

pt ps pu 

+ trm / F"{Z{s)) / Dstrn, / aDiu,r)drdub{s)ds. 
Jo Jo Jo 

Now we collect terms. The two terms on the right-hand side of ()5.6p are the second term 
in (j5.4p and the first term in ()5.5p . Moreover, the last three terms on the right-hand side 
of (|5.7p are equal to the fourth term in (|5.4p , the second term in ()5.5p and the fifth term 
in (|5.4p . Furthermore, the last four terms in ()5.8p are the sixth to ninth term in ()5.4p . 
Finally, the first term one the right-hand side of (|5.8p is twice the last term in (|5.4p . 
Hence, by adding and subtracting the last term in ()5.4p to this equality, we obtain the 
assertions. □ 

Now we provide a sufficient condition under which (j5.2p holds relying on Holder's 
Inequality. 

Remark 5.3. For the term C{s) in l\5.2\) we have to deal with the following three terms 
which constitute its IJ^'P -norm. We will apply Holder's Inequality with respect to oj (note 
that g is assumed to be deterministic) with the two conjugate exponents (?i,(?2 £ [l^oo] 
(where if one of them is infinite, the supremum norm has to be used). We will now 
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assume that Y{s) € -^^2(^2; ^3)1 but a similar calculation would also be possible with 
a{s) e L2{'Hi,Ti2)- We obtain 



E[\\Y{s)g{s,s)a{s)\\l^^^^^^^^]ds 









fT 
/O 

similarly for the Malliavin derivative using ()2.2 



and finally for the second Malliavin derivative 
rT rT rT 



L 



X 

rT 



+ (lE[r(«)Hi%„«3)])'^" [ [ (lE[llAAa(.)|Q,^(^^,^^)^^^,^J)^/'Mtdr)d.. 

From the last equality we see that Y , DY , DDY , a, Da and DDa appear and have to 
satisfy some integrahility conditions with respect to the temporal and spatial argument. 
Therefore we conclude that a sufficient condition for the conditions in (j5.2p is that Y 
and a belong to a Sobolev-type space with respect to all arguments (including the random 
argument) , i.e. in some 'h'^''^ -spaces for q sufficiently large, possibly with a different q for 
each of the two. Moreover, g{s, s) has to be in some -space for p > 4 with respect to 

M[0,T]- 

From here special cases to further manipulate the terms in the equality above are 
thinkable. For instance one can assume that a is independent of the noise (what would 
imply its Malliavin derivatives to be zero), or a might be independent of Y in which 
case the expectations involving Y and a would factor out. An investigation of the term 
A{s) in ()5.2p does not yield to any further conditions, except that also the function 
{t,s) 1-^- %) ^'^^ integrable with respect to X^\[o^t\x[o,t]- 
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A first immediate application of Theorem 15.21 is to derive an Ito formula for X. Note 
that due to ()4.ip Y[s) = id-^j ^-^d therefore 



Xit) 



dX{s)= [ ICg{idn,){t,s)ais)6B{s)+tin, [ {DsJCgiidn,){t, s))ais)ds. 
Jo Jo 

Using Theorem 15.21 we get the following result. 

Corollary 5.4. Suppose that X is defined as in where the conditions in As- 

sumption \3.1\ hold, and that F : ^3 is twice Frechet differentiable. Assume 

moreover that a is twice Malliavin differentiable and that g has the form ()4.5p . Then 
F'{X) € X-^(0,t) for all t G [0,T] and 

F{X{t)) = F{0)+ f F'{X{s))dX{s)-}-tvn, f F"{X{s)){g{s,s)a{s)){g{s,s)a{s))ds. 
Jo ^ Jo 

Note that this is consistent with the Ito Formulas in the one-dimensional case treated 
in [U Theorems 1, 2]. In fact, in that setting we have a = 1, Tii = T-L2 = K, X{t) = 
Jq g{t, s)dW{s) and R{s) := g{s,r)g{s,r)dr. Therefore, Corollary 15.41 implies 

F{X{t)) = F(0) + 1^ (^F'{X{s))g{s, s) + F'{X{u))^{u, s)dv}j 5B{s) 

+ {DsF'{X{s))g{s,s) + j^D,F'{X{u))(t>{u,s)di}jds (5.9) 

F"{X{s))g{s,s)g{s,s)ds. 
Note that DsF'{X{u)) = F" {X{s))g{u, s) for u>s. Note furthermore that 

—^=g{s,s) +2 J g{s,u)(l){s,u)du. 
This and Fubini's Theorem yield that the last three terms in ()5.9p are equal to 

^ 1^ F"{X{s)) (^g{s, sf + 2 g{s, u)4>{s, u)dv^ = I F" {X{s))dR{s). 
This yields 

= F(0) + f F'{Xis))dXis)+ f F"{X{s))dRis), 
Jo Jo 

which is the formula in [U Theorems 1, 2]. 

A second application of Theorem 15.21 is to calculate in the case Ti^ = M. For this 
we suppose the same conditions as in Theorem 15.21 assume Ti^ = R and F{x) = x^. 
Then applying Theorem 15.21 yields 

-{Z{t)f = I Z{s)Y{s)g{s,s)a{s)5B{s) + ( [ Z{s)Y{s)^{s,u)a{u)6B{u)ds 
2 Jo Jo Jo 
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+ trn, [ {DsZ{s))Y{s)g{s,s)a{s)ds + tTn, [' Z{s){DsY{s))g{s,s)a{s)ds 
Jo Jo 

+ trni [ I {DsZ{s))Y{s)(l}{s,u)(j{u)duds 



10 Jo 

tr-Hi I I Z{s){DsY{s))(l3{s,u)a{u)duds 
Jo Jo 



~ll \\yis)g{s,s)a{s)\\l^^^^^^^^ds. (5.10) 

If we now introduce the symbolic notation 

dZ{s) =ICg{Y){; s)a{s)6B{s) + tr^, i?s/Cg(y)(-, s)a{s)ds 

=Y{s)g{s,s)a{s)6B{s) + [ Y{s)^{s,u)a{u)6B{u)ds 

Jo 

+ trm {DsY{s))g{s, s)a{s)ds + tr^, / {DsY{s))^{s, u)a{u)duds, 

Jo 

then we see that the term Z{s)dZ{s) accounts for the first, second, fourth, fifth and 
sixth term (use the rule in Proposition I2.2p in ()5.10p . The third and last term remain 
as correction terms and one has the formal result 



t ft 

Z{s)dZ{s) + tvu, / {DsZ{s))Y{s)g{s,s)a{s)ds 
Jo 



\Y{s)g{s,s)ais)\\%^^^^^^)ds. (5.11) 



One could also alter the correction terms. For instance, one could also expand the second 
term on the right-hand side of (j5.1ip by calculating DsZ(s). This would reverse the sign 
of the third term on the right-hand side of (|5.1ip . but one would get other correction 
terms. 



6 A random-field approach to the X-integral 

In this section we investigate whether the random-field approach pioneered in [19] gives a 
reasonable interpretation of the X-integral, i.e. whether we can derive an integral which 
has the form f^^ y{s, y)X{ds, dy), where 

X{t,x)= f I g{t,s;x,y)a{s,y)M{6s,dy), (6.1) 
Jo Jw 

where, similar to the terms in (jl.ip and (|1.2p . g is a deterministic function, a and Y are 
random fields and M is a martingale measure and the integral in (j6.ip is understood in 
the Walsh sense. In the first part of this section we quickly review the concept of Walsh 
integration and summarize the main ideas of a minor generalization for anticipating 
integrands. In the second part we present the random-field X-integral in the special 
case of homogeneous noise and show that this integral and the one derived in Section [3] 
coincide. 



31 



6.1 Walsh integration 



Let (Mt{A);t G [0,T],A € ^^(M )) be a worthy Gaussian martingale measure, where 
are the bounded Borel sets. This means that each Mt{A) has a Gaussian distribution 
on M. Fix the filtration {J't)t&[o,T] ^ the one generated by the martingale measure and 
augmented by the sets of probability measure zero. The martingale measure takes values 
in L^(r2). Since the martingale measure is a martingale in t for all fixed sets A G Bi,{^'^), 
we can associate a quadratic covariance functional to it, denoted by QMit, A, B) = 
{M{A), M{B))t. We can define a stochastic integral of elementary processes f{s,x,u;) = 
^{a,b]{s)lA{x) Z (lo) , where 0<a<b<T,AG Bi){M.'^) and Z is a bounded Ja-measurable 
random variable, with respect to the martingale measure by 



f{s,y)M{ds,dy) = {Mt^b{A) - Mt^a{A))Z, 

Jr'' 

and for simple processes (linear combinations of elementary processes) by an obvious 
linear combination. Then, as in Ito integration, this stochastic integral is extended to a 
larger class of integrands by using the isometry 



E 



f{s,y)M{ds,dy) 



E 



JR'^ JR'^ 



f{s, x)f{s, y)QM{ds, dx, dy) 



(6.2) 

where we extended the covariance functional to a measure on [0, T] x M'^ x W^. In order to 
be able to do that, one needs the fact that the martingale measure is assumed to be wor- 
thy. This implies that there exists a dominating measure Km-, that is a positive definite 
measure with some regularity conditions for which |Qa./([0, t], ^, -B)| < KM{{Q-,t\, A,B) 
for all t G [0,T] and A,Bg Bi^iW^). In this case one gets an upper bound for the 
isometry (|6.2|) 



E 



/(s, x)/(s, y)QM{ds, dx, dy) 



< E 



JR'' 
2 



\f{s,x)\\f{s,y)\KM{ds,dx,dy) 



The norm || • ||o is actually induced by a scalar product and one defines two possible 
domains and Vq of the stochastic integral. The former one is defined to be the set 
of all functions / with ||/||+ < oo and the latter one is defined as the completion of the 
space of simple processes with respect to the scalar product (•,-)o- Note that by [191 
Proposition 2.3], is complete and the simple processes form a dense subset in this 
space. Vq on the other hand is a Hilbert space, but does not coincide with the set of 
all processes such that ||/||o < oo. Note that for particular choices of M, there might 
be even Schwartz distributions in T'o- A particular example for this construction is the 
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case of spatially homogeneous noise, where 



QM([0,t],A5) =i^M([0,t],Ai?) = r [ [ lio,t]{s)lA{y)lB{y - x)dyT{dx)ds, 

Jo JRd JRd 

(6.3) 

for a nonnegative, nonnegative definite tempered measure T. Then one can find a spectral 
representation of the || • ||o-norm so that in this case Vq contains distributions. Note 
however, that in this article we only use the integration concept by Walsh for functions 
since we are exclusively interested in the case when g' is a function, as this is the case 
with ambit fields in Example ll.il However, in Section [7] we will also admit distributions 
for g which calls for the Dalang integral [B]. 

For a correct treatment, we need to provide a certain extension of the Walsh integral 
in the case of homogeneous noise so that it can handle anticipating integrands. This is 
however not very difficult if one keeps in mind the reformulation of the Walsh integral as 
a sum of independent Brownian motions on a Hilbert space, see [5] or as the divergence 
operator of Malliavin calculus, see [18] . 

We define Ti to be the completion of the simple functions on by the scalar product 

{f,9)H= [ U * g){z)T{dz) = f [ fiy)giy-z)dyridz), 

with r as above and T-Lt ■= L'^{[0,T];'H). In contrast to [6] we do not need to introduce 
the spectral representation of this scalar product since we need / and g to be functions. 
Then the worthy martingale measure M can be extended into an isonormal Gaussian 
process {M{h);h G T^j') with respect to which we can do Malliavin Calculus. With 
this we can define B'^'^{n) and the important space h^'^{n) := L^{[0,T];O^'^{n)), with 
norm 

On this space the divergence operator 6 : L?'{Q.; T-Lt) L?'{^) is well-defined and contin- 
uous, see [m (1.47)]. For simple processes /, which are all in L^'^('H) the anticipating 
Walsh integral is given by 

d d 

/ fit, x)M{St, dx) = J2 Zj {MtAb, (Aj) - Mt^a, (Aj)) + Y.{DZ„ l[,^,6^.]U,)w^. 
j=i j=i 

Then one can extend this integral to integrands in L^'^(?^), where the Malliavin dif- 
ferentiable simple processes are dense in. This extension in done using the following 
isometry 



E[(n-M)^ 

////// Ds2,y2-Z2u{si,yi) 

Jo Jo JW JR'i JW JR'i 



E[\\ufn^]+m\Dn\\ 



m\u\\i^]+^ 
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X Ds^^y^^z-^u{s2,y2)dyidy2T{dzi)T{dz2)dsids2 



(6.4) 



Another way to think about this integral is by using the equivalence of the anticipating 
Walsh integral and a sum of anticipating integrals with respect to Brownian motion. In 
fact we have for all u G h^''^(T-L) 



where {ek)k£N is a CONS of T-L, B is a Brownian motion on H and B'' := {B, e^)-^ are in- 
dependent real-valued Brownian motions, for which anticipating calculus is well-known. 
Therefore, one could also take the last equality as the definition of the anticipating Walsh 
integral. 

6.2 The random-field X-integral 

Now we start with the definition of the X-integral using a random-field approach. In 
contrast to the definition of the X-integral in Section [3] we have to start here with 
defining the stochastic integral on elementary processes first and then extend it to simple 
processes and further. This is due to the fact that a similar integration by parts procedure 
done in Section 3 is not easily applicable since one would get nontrivial boundary terms 
which are hard to interpret. 

Take H and T-Lt to be the Hilbert spaces defined in the previous section. The mar- 
tingale measure that we will use in the following will be the one using the homogeneous 
noise in ()6.3p . which appeared in [6|. We will however not use it in its full generality, 
but we only look at functions / for which the following norm 



This is a norm - and in fact it is generated by a scalar product (•, •)o - as one can see 
by going to its spectral representation. We will however not do this, since we can not 
yet deal with distributions as integrators g in the kernel p.4p . Therefore, we define the 
Hilbert space Ht as the completion of the simple functions with respect to the (•, •)o- 
This is then a function space which does not contain distributions. We also look at 
random functions in L^(0;'Ht) equipped with the obvious norm and scalar product. 
The following assumptions are made: 

Assumption 6.1. Fix T > and let t e [0,r], and let M = {Mt{A);t G [0,r],y4 G 
B},(^'^)) be a worthy martingale measure with covariation measure (16. 5p . Furthermore, 
let g : [0,T]^ x (M'^)^ — > M 6e a deterministic function for all < s < t < T and 
all x,y £ and let {a{t,x);{t,x) G [0,r] x M'^) be random field such that {s,y) i-^ 
g{t,s,x,y)a{s,y) is integrable with respect to M for all {t,x) G [0,T] x M'', i.e. that this 
random function is in L^'^(?^). Assume for all (s,y) G [0,T] x that the function 
{t,x) I— 7- g{t,s,x,y) has bounded variation on [u,v] x M'^ for all < s < u < v < t. 



[ [ u{s,y)M{5s,dy) = y, I {u{s,-),ek)H5B\s) 
Jo Jr" f^^^Jo 




(6.5) 
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Now we turn to the definition of the integral. To this end, let first 0<a<6<r, yl€ 
be of the form A = X'^^i[aj,bj] and set Y{t,x,u;) = l(^a,b]{'l^)^A{x)Z{u), where Z 
is a bounded Malliavin differentiable random variable, not necessarily measurable with 
respect to Ta- Note that since the spatial argument of the integrator process X only 
appears in g, and since X is linear in the kernel g, one can derive using ()2.ip 



X{t,A)= [ [ g{t,s,A,y)a{s,y)M{ds,dy) 

Jo JRd 

^A{z)g{t, s, dz, y)a{s, y)M{ds, dy). (6.6) 





In the calculation that follows we aim at arriving at a similar kernel as in ()3.4p . In order 
to keep the notation tidy, we do this calculation for the case when g{s,s; •, •) is a well 
defined object, thus yielding a kernel as in (|3.5p . By following the same arguments as 
given below, one can however also arrive at ()3.4p . Using the obvious definition for the 
stochastic integral for simple processes and (j6.6p . one derives 

f [ Y{s,y)X{ds,dy) 
Jo JR-i 

:= Z{X{t ^h,A)- X{t A a. A)) 
ft 



Z J^^ ( l^^lAiz)g{s,s;dz,y) 



+ l[o,b](s)( / 'i^Aiz)g{b,s;dz,y) - lAiz)g{s,s;dz,y) 

\ JRd JR-i / 

-\oa{s)\ / lA{z)gia,s]dz,y) - lAiz)g{s, s;dz,y)] \a{s,y)M{6s,dy) 

KJRd JRd J J 

^1 I \\a,b]is) I lA{z)g{s,s;dz,y) + liQ^a]{s) [ [ lA{z)g{du,s;dz,y) 



+ 1m(s)/ / '^Aiz)gidu,s;dz,y)]a{s,y)M{6s,dy). 

Js JR-i J 

Note that the last two integral terms one the right-hand side of the previous equation 
are equal to 

f'bAt r r 

I I lAiz)gidu,s;dz,y) = / l[a^f,]{u)lA{z)g{du,s;dz,y). 

JaVs JR'' Js JR''- 

Plugging this into the equation above yields 

/* / Y{s,y)X{ds,dy) 
Jo Jr'' 

l[a,b] {s)lA{z)g{s, s; dz, v) + Ma,b] {u)lA{z)g{du, s; dz, y)^ a{s, y)M{6s, dy). 



j-t 

Z 



(6.7) 
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Similar to (j3.4p . we define an integration kernel 
)Cg{h){t,s,y) := h{s,z)g{t,s;dz,y) + / {h{u,z) - h{s,z))g{du,s;dz,y). (6.8) 



Note that this term is well-defined if is a function of bounded variation spatial argu- 
ment z and of bounded variation on subintervals bounded away from s in the temporal 
argument t. As in Section[3l one can rewrite this kernel in some situations. If g{s, s; z, y) 
exists almost everywhere (with respect to dyT{dz)ds), then one has 



ICg 



{h){t,s,y)= h{s,z)g{s,s;dz,y) + / h{u, z)g{du, s;dz,y). 

JR'' Js 



If the function g has {d + 1) partial derivatives, then 

f d'^ 
K:g{h){t,s,y) = / h{s,z)- —g{t,s;z,y)dz 

Jjjd OZi . . . OZd 

r-t r- Qd+l 

+ Is L " 5.9.1... 9./^"' 

Using the kernel in (j6.8p . we obtain by pulling the random variable inside the stochastic 
integral 

* / Y{s,y)X{ds,dy) 

Z [ [ )Cg{l[a,b]U){t,s,y)a{s,y)M{Ss,dy) 

Jo JR'i 

/ K^gi^a,b]'^AZ)it,s,y)a{s,y)M{6s,dy) + {DZ,}Cg{l[a,b]Wii^ •> *)o-(-, 
JR'i 

The second term of the previous equality can be manipulated as follows in order to get 
an expression which explicitly involves the original integrand Y 

(z)z,/Cg(i[,,6]U)(t,-,*M-,*))wT 

r-t r r 

Ds,y-zZJCg{l[a,b]'^A){t, s, y)a{s, y)dyT{dz)ds 
III Ds,y-,]Cg{Y){t,s,y)a{s,y)dyT{dz)ds. (6.9) 



JR"^ 
t 



In total, we have for the random-field X-integral 

' Y{s,y)X{dsAy)= f I }Cg{YKt,s,y)a{.s,y)M{6s,dy) (6.10) 



Jr'^ Jo 



+ f I I D,^y^,lCg{Y){t,s,y)a{s,y)dyT{dz)ds. 
Jo JR'^ JR'' 
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First we see the similarities to the definition in (|3.6p . where here we add the spatial 
component which in turn forces us to include a possible spatial correlation into the 
definition of the integral. Note that if we consider the uncorrelated case, then F = 5o 
and one gets an expression without z. Now we extend the definition of the integral by 
using the isometry in (j6.4p so that the integral (j6.10p is defined for all Y € L^'^(?^) for 
which similar integrability conditions as in Definition 13.21 hold. 

Definition 6.2. Let X be the random field defined in ()6.ip together with Assumption \6. 1[ 
We say that the random field Y = (Y{t, x); [t, x) G [0, T] xM'^) € Vq belongs to the domain 
of the stochastic integral with respect to the random field X, if 

1. the process {Y{s, z))^^T^d is integrable with respect to g{t, s; dz, y) almost surely and 
{t, s,y)-almost everywhere, 

2. the process 2) — ^(s, 2:))„g(5 jjxRd is integrable with respect to g{du, s;dz,y) 
almost surely and {s,y)-almost everywhere, 

3. {s,y) I— )• }Cg{Y)(t, s,y)a{s,y)l[Qij{s) is in the domain of the martingale measure 
M, i.e. /Cg(F)(t,-,*)c7(-,*)l[o,t](-) G and 

4. lCg{Y){t, s,y) is Malliavin differentiable with respect to Dg^y^z for all s G [0, t] 
and y,z G and the random field {s,y) i-> t^i-u^ Ds^y-z{^g{Y){t, s,y))a{s,y) is 
Lebesgue integrable on [0, t] x almost surely. 

We denote this byY^ 2^'''" ([0, t] x W^) and the integral Jg J^^ Y{s, y)X{ds, dy) is defined 
by dSIOD. 

The main reason why we have introduced this integral is to be able to integrate 
processes {Y{t,x,y);t G [0,T],x,y G M'*) with respect to t and y and obtain a random 
field which has a pointwise (in x) interpretation as a random variable, rather than as an 
element in some abstract Hilbert space. This could serve to deduce properties of ambit 
fields such as path continuity, or existence of densities at each point. 

Next we could easily derive similar properties for the random-field X-integral as in 
Section 14. H but we restrict ourselves to showing that the random field integral can be 
rewritten as a Hilbert- valued stochastic integral when we interpret the process X as a 

valued Wiener process, as it was shown in [51 [T5]. 

Proposition 6.3. LetY e 1^ {[0,t]xR'^). ThenYeI^{0,t) with Hi = V. and = R 
and ^ ^ 

f f Y{s,y)X{ds,dy)= f Y{s)dX{s), 
Jo JR-i Jo 

where X is interpreted as a random field on the left-hand side and as an 'H2-valued 
Volterra process on the right-hand one. 

Proof At first we argue that if F G X^([0,t] x M''), then Y G l^{0,t). In fact, from 
the remark at the end of the previous section it follows that if Y is in the domain of the 
martingale measure integral, then it is also in the domain of the integral with respect 
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to a ?^-valued Brownian motion. Furthermore, the properties I6.2IBI H1 are exactly 13. 211 13] 
written out in the special case of Tii = %. This yields the first assertion. 

Secondly we show that the stochastic integrals with respect to the martingale measure 
and the 7^-valued Brownian motion are equivalent. This follows by a straightforward 
adaption of the proof of [5', Proposition 2.6], because by Assumption 16 .251 we conclude 
that 

lCg{Y){t,;*)a{;*)\^^t^{-) e L\^-Ut) 

and 



t 

JR'i 



lCg{Y){t, s, y)a{s, y)M{5s, dy) = Y. f ^' *)^(^' *)' ek)^dB\s), 



where B is an ?^-valued Wiener process and := {B,ek)n- This yields the equality of 
the stochastic integrals in the X-integral and the random-field X-integral. 

For the pathwise integrals we start with the elementary processes given hyY(t,x,u}) = 
^[a,b]{.'t)^A{x)Z{uj) and use the representation in (|6.9|) to obtain 



t 

JR'' 







Ds,y-zK,g{Y){t, s, y)cr{s, y)dyT{dz)ds 
ft 

{DsZ, }Cg{lia,b]'^A){t, s, *)a{s, *))nds 
= X] / i^sZ,ek)n{}<^g{^a,b]'^A){t,s,*)a{s,*),ek)Hds 

= XI / {Ds,klCg{Y){t,s,*)a{s,*),ek)n<is 

= tr^ / DsJCg{Y){t,s)a{s)ds. 
Jo 

Then we extend the last equality to all processes Y £ Vq which satisfy Definition 16.2111 
This, together with the equality above, implies the assertion. □ 

Now we treat a first example for the random field X-integral, another one will follow 
in Section [71 Here we focus on the nonlinear stochastic heat equation with null initial 
condition 

d \ 

— (t,x) - A ]u{t,x) = b{u{t,x)) + a{u{t,x))F{t,x) {t,x G (0,r] x M.'^) 



dt 



u{0, x) = X € 



nd 



The fundamental solution of the associated PDE (the heat equation) is given in any 
spatial dimensions by 

, , 1 r-\x\^' 

9{t, X) = ,^ exp 
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Let in the following without much loss of generality d = 1. Note that the mild solution 
to the SPDE above is given by 



u{t,x)= / g{t-s,x-y)a{u{s,y))M{ds,dy)+ / g{t - s, x - y)b{u{s,y))dyds, 
Jo JR'' Jo Jr"^ 

with M being the martingale measure corresponding to the noise F. We see that the 
stochastic integral has the form of a Volterra process and therefore we can apply the 
integration theory developed in this section. Since the fundamental solution g is differen- 
tiable in both arguments and stationary in time and space, we can calculate }Cg{Y)(t, s, y) 
to be 

lCg{Y){t,s,y)= [ Y{s,y + z)g{t-s,dz)+ [ [ Y{u + s, z + y)g{du,dz) 

jRd Jo jRd 

= Y{s,y + z)dzgit - s,z)dz + / Y{u + s, z + y)d? ^^g{u, z)dzdu, 

jRd Jo J^d 

where dz denotes the partial derivative with respect to the spatial argument, and d? n 
is the mixed spatial and temporal argument. Note the heat kernel estimates for the 
derivatives 

\d,dtg{t, s, x,y)\< c{x -y){t- ,)-(3+2a)/2 exp ( - ^i^^ 



for a G {0, 1}. This implies that all Y which are almost surely polynomially bounded are 
integrable with respect to g. This implies [1] and [2] in Definition 16.21 In order for [3] to be 
satisfied one needs that lCg{Y){t, ■, *)cr(-, *) € Vo- This depends on the concrete form of 
a and F. In order for the pathwise integral to well-defined in L'^{Q) which implies IU one 
has to assume that Y is Malliavin differentiable and that the integrand of the pathwise 
integral is in L'^{[0,T] x n;L^{R'^)). 



7 Generalization of the random-field X-integral 

In this section we follow two objectives. First, we want to present another explicit exam- 
ple for the X-integral defined in Section [3] and thus show how nonregular fundamental 
solutions to partial differential equations enter into this framework. Secondly, we want 
to show a first idea how to generalize the random-field X-integral in Section [6] in a way 
similar to the generalizations of the Walsh integral by Dalang in [6] , going from functions 
to distributions as integrands. More specifically, we investigate under which conditions 
on Y the kernel ICg{Y){t, s) is well defined if g is chosen to be the fundamental solution 
of a partial differential operator which is singular in the spatial argument. 

For this we focus on the specific example of a wave equation in different spatial di- 
mensions d given by 

(^^-Ad^u{t,x) = 6o, (7.1) 
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with null initial conditions, where A^^ is the d-dimensional Laplace operator and 60 is the 
Dirac distribution. The solution to this equation is the the fundamental solution to the 
wave equation, which differs with the spatial dimension. In the dimensions d = 1 and 
d = 2 the fundamental solution is a function, which we can treat with similar methods 
as in the example with the heat kernel in Section [6l 

For d = 3 the fundamental solution is given by g{t, s) = cp^_g/ {t — s), where c > is 
a constant and is the surface measure on the sphere in three dimensions with radius 
t. Therefore, the mass of this measure is equal to ct^. Since this fundamental solution 
is not a regular function anymore, the random-field X-integral of Section [6] is no longer 
well-defined. So we first use the Hilbert- valued integral from Section [3] with the Hilbert 
space T-L introduced in Section 16.11 Let y be a random element with values in the linear 
functionals on Ti. In order to answer the question about the form of fCgiY) we perform 
integration by parts to get 

-i/C,(y)(t,s) = Y{s)^^ + {Y{v + s)- A^dt; 

Y'{v)^dv 

V 



3 3 

= Y{s)^ + iY{v + s)- Y{s))^ 

t — S V 
„3 pt-s 3 

= y(t)i^ _ / Y'{v)(^dv, 
t- s Jo V 

where the integral in the last line is understood as a Bochner integral. Since the total 
mass of the measure grows in v as v'^, the last equality tells us that }Cg(Y)(t, s) is 
well defined as an operator if Y is differentiable and does not have a singularity of order 
greater or equal than 2 at zero, i.e. Y € o{v~^) at zero. In this case, ICg{Y){t, s) exists 
as a linear operator on T-L defined by 

ICg{Y){t, s)f{s) = _ Y'{v)^^dv, 

t — S Jq V 

where 

Plifis)) = [ f{s,y)plm = v' [ f{s,vy)pl{dy). (7.2) 

Under these conditions the X-integral is given by 

rt rt / 3 rt — s 3 \ 

j Y{s)dX{s)=j [Y{t)^^-J Y'{v)^dvy{s)dB{s) 

+ trn 1^ [DsY{t)^^ - 1^ ' DsY'{v)^dv^ a{s)ds, 

and this holds for any a for which (|7.2|) is finite and for which the integrability conditions 
from Definition 13.21 hold. 

The same procedure can be done in even higher dimensions d. The fundamental 
solution to the wave equation becomes in these cases 



git) 



r{d/2)^t>o[tm) T ltd IS Odd 

lT(dm^>°^tdi) i^ -m)+ ifdiseven 
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This can also be treated by integration by parts. Let us show this for the case d = 4, 
where all the higher-order cases follow in the same manner. For d = 4 we deal with the 
regular distribution (t^ — | • given for each test function / by 

JB4o,t) V* - \y\ JBi{o,i) ^/l - |yr 

Hence, 

c^^)Cg{Y){t,s) 

= Yis)git -s) + fj^ [Y{v + s)- Y{s)) ^ Q (,2 _ | . |2)-i/2^^ 

I — S Jo ^ 

(7.3) 



In order for this to well defined we need the following condition for an integration by 
parts procedure 

limI^(.2-|.|V/^ = 0. 

This implies that Y must be differentiable in the temporal argument and Y G o{v~^) at 
zero, since the expression in the brackets is of v'^. In conclusion, in order for the integral 
in ()7.3p to be well defined, we need Y to be twice differentiable and Y' € o{v~'^) and 
Y" G o(u~^) at zero. 

Now we want to go a step further and show that even the random-field integral exists. 
We show this in the case d = 3, but the method can be applied in higher dimensions, 
too. So we assume in the following that Y{s,*) vanishes at infinity for all s G [0,r] 
and that it is differentiable in s and in x. Then we can apply an integration by parts 
procedure and obtain 



Y{s,z)g{t,s;dz,y) = / g{t,s;z,y)d^Y{s,z)dz, 

where = d"^ /dzidz2dz3. Under these conditions we can rewrite the integral with 
respect to the martingale measure and use its isometry to obtain 



E 



ICg{Y){t,s,y)a{s,y)M{ds,dy) 
JR'' 



E 



t 

JR'^ 



ICg{Y){t, s,y - z)ICg{Y){t, s, y)a{s, y - z)a{s, y)dyT{dz)ds 

Note that by integration by parts with respect to the temporal and spatial argument 
ICg{Y){t,s,y-z)}Cg{Y){t,s,y) 



(7.4) 
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^ p pt—s p 

/ dliY{t,x^)g{t-s;x^-y + z)dx^+ / 8^ ^,Y{v,x^)g{v-dx^ - y + z)dx^&, 

Jr3 Jo Jr3 ' 

( [ dl2Y{t,x^)g{t-s;x^ -y)dx^+ [ [ ^,Y{v,x^)g{v;x^ - y)dx^dv] . 
\Jr3 Jo Jr3 ' J 



(7.5) 

Factoring out this product and substituting this into (|7.4p yields to four stochastic 
integrals of which we will explicitly only treat the first one. So we have 



E 



/ / / [ dliY{t,x^)g{t-s,x^ -y + z)dx^ [ dl2Y{t,x'^)g{t-s,x'^ -y)dx'^ 

Jo J«.3 Jk3 Jk3 Jir3 

X cr{s,y — z)a{s,y)dyr{dz)ds 
< [ [ sup E[iey(t,xi)53,y(t,x2)||a(r,y)|2] 

JR3 JR3 {r,j/)G[0,T]xR3 
t r 

exp{i^{x^ + x^))\\^A{t - s){0\'^fi{d^)dsdx^dx^ 



10 JM? 

<[ {E[{d'^Y{t,x)fP])^^^dx sup E[cj(r,y)2'?]i/^ / / |^A(t - s)(OlV(de)ds 

where in the last line we used Holder's inequality with p^^ + q^^ = 1. The other three 
terms in (|7.5|) and the pathwise integral involving the Malliavin derivative can be treated 
in a similar way, so that at the end the random field integral f^3Y{s,y)X{ds,dy) 
exists when g is the fundamental solution of the three-dimensional wave equation if a 
has uniform (in time and space) 2g-moments, if 

/*/ |^A(t-s)(0|V(d0d5<oo, 
Jo Jr^ 

for being the Fourier transform, and if 

{E[idMt,x)f^]f^''^ + [ {'E[idl^Yiu,x)f^]f'dxdu < oo. 

Jo Jr^ 

Here we see a certain trade-off between the integrability conditions of Y and a, if we 
let one of the two to be more irregular, the other one has to be more regular. A similar 
condition with the Malliavin derivatives of Y has to hold so that the pathwise integral 
is well defined. All this means that Y has to be in some Sobolev-type space where the 
(weak) derivatives of the process and of the Malliavin derivatives of the process exist 
and an integrability condition has to hold. These conditions are far from optimal, but 
serve as a first step towards a general existence theorem for the random-field X-integral 
when g is a distribution. 
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